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I. INTRODUCTION 

After decades of theoretical efforts, understanding the 
low-energy properties of bidimensionnal frustrated Quantum 
Heisenberg Antiferromagnets (QHAF) is still a notoriously 
puzzling problem. In contrast to unfrustated models where 
powerful paradigms have emerged to describe both the mag- 
netically ordered ground state (Neel state) and the excitations 
(magnons), a general theoretical framework is still lacking. 
Not only the low-energy physics is in general poorly char- 
acterized - one can have in mind the lack of consensus on 
archetypal models such as Ji — J2 or kagome antiferromag- 
nets - but generic tools to identify the low energy degrees of 
freedom and understand their effective interaction are miss- 
ing 1 . Recent discoveries of new materials such as herbert- 
smithite 2 , believed to be almost perfect spin- 1/2 kagome an- 
tiferromagnets, have revived the former interest for bidimen- 
sionnal frustrated QHAF and call more than ever for novel 
theoretical approaches that could provide a quantitative under- 
standing of these materials and their exotic low-temperature 
non-magnetic phases 1 . 

The underlying hardness of the problem comes from the 
fact that the usual theoretical frame for treating such systems 
typically requires to (i) identify a well controlled limit (clas- 
sical limit, unperturbed limit, . . . ) and then (ii) introduce 
quantum fluctuations and/or corrections to this limit in a pu- 
tatively perturbative way. This scheme however fails for bidi- 
mensionnal frustrated QHAF since none of the above men- 
tioned points is fulfilled : (i) because of frustration the classi- 
cal limit is highly degenerate and often not fully understood, 
(ii) quantum fluctuations generally do not act as a small per- 
turbation on the classical limit. As a direct consequence, the 
low-energy physics of these systems is characterized by very 
exotic magnetically disordered states such as Resonating Va- 
lence Bond (RVB) states 3 , Valence Bond Crystals (VBC) 4 or 
spin liquids 5 . 

The aim of this article is to build a non-perturbative scheme 
allowing to derive an effective low-energy Generalized Quan- 
tum Dimer Model (GQDM) that captures the physics of the 
microscopic model. In turn, the GQDM can be used to ef- 
ficiently investigate the low-energy properties of the micro- 



scopic Hamiltonian either by numerical techniques on large 
systems (by Exact Diagonalizations or Quantum Monte Carlo 
when the GQDM has no sign problem) or by analytical tech- 
niques such as gauge theory mapping. 

Basically, the method relies on (i) the identification of a ver- 
satile and relevant manifold of states at low-energy and (ii) the 
projection of the microscopic Hamiltonian in this manifold. 
Such a framework, first initiated by Rokhsar and Kivelson 6 in 
a different context, has been put here on a formal basis and 
greatly extended, providing a novel versatile and systematic 
expansion scheme to deal with frustrated magnets. 

The paper is organized as follows: in the first part we 
present the derivation scheme focusing on important results 
and highlighting practical implementation rules. To increase 
readability, we make an extensive use of appendices in which 
all technical details and demonstrations are postponed. The 
second part is devoted to the application of the method to 
the kagome antiferromagnet. We derive the parameter-free 
GQDM Hamiltonian that has been proved to give deeper in- 
sights to the low energy physics of the kagome antiferromag- 
net in a companion paper 7 . Then we discuss its key properties 
focusing on the fact that this untuned GQDM model lies at 
the vicinity of several competing exotic phases, a Z2 dimer 
liquid and two VBC phases, sheding light on the critical prop- 
erties of the microscopic model. For completeness, we pro- 
vide complementary results for the Kagome lattice at [URL 
will be inserted by AIP]. 



II. EFFECTIVE HAMILTONIAN MAPPING SCHEME 

A. Low energy manifold 

In the past decades numerical studies played a central role 
in establishing unbiased and reliable results on bidimension- 
nal frustrated QHAF 1 . More precisely, Exact Diagonaliza- 
tions (ED) of the Heisenberg Hamiltonian on finite and gener- 
ally small clusters were able to provide precise insights on 
these problems. Other numerical methods, such as Quan- 
tum Monte Carlo (QMC) or Density Matrix Renormalization 
Group (DMRG), have proven their relative inadequacy to treat 
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such systems and failed to be generalized for reasons that ap- 
pear more and more to be fundamental rather than purely tech- 
nical. 

Based on ED results, two salient facts can be retained to 
draw a phenomenological low-energy picture: bidimension- 
nal frustrated QHAF (i) remain magnetically disordered at 
zero temperature (i.e. the zero temperature spin-spin corre- 
lation function (Sj.Sj) does not display long range antifer- 
romagnetic order), (ii) do not break the SU(2) symmetry, the 
groundstate and low-lying excitations are singlets (i.e. eigen- 
states of the total spin S with a zero eigenvalue). 

Singlet states are intimately related to Valence Bond (VB) 
states, namely products over pairs of spins of arbitrary range 
SU(2) singlet wavefunctions (also called SU(2) dimers in the 
following). Indeed, any singlet state can be expressed as a 
linear superposition of VB states 8 ' 9 . It is therefore tempting 
to build a low-energy singlet theory by representing the orig- 
inal spin Hamiltonian in the VB basis rather than using spin 
variables. However, the set of all VB states is massively over- 
complete, as the number of VB states is much larger than the 
number of singlets. This originates from the fact that in the 
VB representation dimers are allowed to connect arbitrarily 
distant sites. On the other hand, Liang, Douc,ot and Anderson 
showed that a deep connection exists between the spin cor- 
relation length of a singlet wavefunction and its bond-length 
distribution 10 : a finite correlation length for (S,.Sj) can only 
be obtained for a sufficiently fast decay with dimer length of 
dimer amplitudes in the wavefunction. 

This strongly suggests that a versatile framework for de- 
scribing the magnetically disordered low energy manifold of 
bidimensionnal frustrated QHAF can be obtained by consid- 
ering only the restriction of the VB states to short range or 
Nearest Neighbor Valence Bond (NNVB) states. The pre- 
cise numerical check of this statement as well as the ability 
of the NNVB set of states to capture the low energy physics 
of bidimensionnal frustrated QHAF has been performed in 
previous work in the case of the J\ — — J3 model on a 
square lattice 11 and for the kagome antiferromagnet 1213 . At 
this point, it is important to emphasize that such a direct eval- 
uation of the quality of this approximation is an essential step 
of the method. Indeed, the projection approximation and the 
analytical scheme presented in the following sections are dis- 
tinct issues as the latter does not provide built-in control on 
the validity of the former. 

As a concluding remark let us mention that, contrary to the 
unrestricted VB set of states, the NNVB states are linearly in- 
dependent for most low-connectivity lattices. This can be nu- 
merically checked for reasonably large systems on the square, 
triangular and kagome lattices 14 and has been recently proved 
for the kagome lattice 15 . However, as the number of NNRVB 
states on those lattices is smaller than the total number of sin- 
glets, it is clear that this restriction further reduces the singlet 
Hilbert space which is a key advantage for numerical compu- 
tations. 



B. Effective Hamiltonian 

For convenience, most of the illustrations provided in this 
section are represented using the kagome lattice. However, 
the formalism is general and all technical details, provided in 
the Appendices V A to V H and to which we extensively refer, 
use a lattice-independent formulation. 

1. Linear algebra problem 

Non-orthogonality . A crucial property of VB states is their 
non-orthogonality: for any \<p) and the scalar product 
(<p\ip) is always non-zero and, as recalled in Appendix V A, 

=a N - 2ni( -^\ (1) 

where N is the number of sites of the system and ni (ip, ip) 
is the number of loops in the overlap graph (see figure 1). 
In the general case, the overlap (y>|^>) is a signed quantity. 
This is a direct consequence of the SU(2) dimer wavefunction 
antisymmetry which requires in turn to specify a conventional 
orientation of lattice bonds. 

When restricting to NNVB states, this choice is constrained 
by the nature of the lattice. The bosonic convention requires 
a prescription of bond orientations and hence is more adapted 
to bipartite lattices. On the other hand, the fermionic con- 
vention, being free of any dimer ordering, is convenient for 
non-bipartite lattices (see Appendix V A). Using the fermionic 
convention (for arbitrary lattices) or the bosonic convention on 
bipartite lattices allow to absorb the sign of (<p\ip) in the single 
parameter a by taking otb = 1/ \pl (respectively a/ = i/y/2). 




FIG. 1: (Color online) Overlap (ip\ip) between two VB wave func- 
tions \ip) and \ip). The amplitude of (<p\ip) is driven by the number 
of loops appearing in the overlap graph (bottom frame). 

Dual basis. When working with a non-orthogonal basis 
{ I ip)} as the NNVB states, it is generally useful to introduce 



its dual basis {|"0*)}> such that 

(ip*\ip) = S^. 



(2) 



Of course, both span the same space, but the dual states are 
clearly not VB states, as the overlap between VB states is al- 
ways non-zero. As we want to work with the NNVB basis, 
one can then further introduce the operator 



6 = J2\x)(xl 



(3) 



that transforms the dual states into NNVB states: 0\ip*) — 
\ip). This operator always exists and it maps any state to a 
state in the NNVB subspace. If the NNVB states are linearly 
independent, which is equivalent to saying that they form a 
basis or that the determinant of O does not vanish, then O is 
invertible within the NNVB subspace. Its inverse is given by 



o- 1 = J2\x*)ix*\ 



(4) 



and transforms the NNVB states back to the dual states, and 



thus the dual states are only defined if O 
it follows, that {^O' 1 ^) = {<p*\i>) = 



exists. In this case 



and {\ip x )} 



{0-V 2 |^} is easily seen to be an orthonormal basis 6 . 

In order to calculate matrix elements in a non-orthogonal 
basis, one may define 



A 



(<P*\A\i>) 

(via- 1 Aw). 



Setting n = 6n H , we find 

H v ,ip = (ip\Hn\ip). 



(5) 



(6) 
(7) 



O v is called overlap matrix and the Heisenberg Hamilto- 
nian for bidimensionnal frustrated QHAF in its general form 
is given by 



(8) 



Matrix elements. The matrix elements can be ob- 

tained from the overlap matrix elements straightfor- 
wardly, generalizing Ref. 16. It can be shown that, 



(9) 



where efy can be derived from the overlap graph. The details 
are provided in the Appendix V A, and following the notations 
of Figure 2, three cases occur: 

a. i and j are two sites lying at an odd distance on the same 
loop in the overlap graph of \<p) and |?/>) and efj = 
-3/4, 
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FIG. 2: (Color online) Application of Pi = 2S,.Sj + 1/2 on a 
loop diagram (<p\ip) for the three types of bonds / = (dashed 
lines) described in the text: "internal odd" (a), "external" (b) and 
"internal even" (c). The state \ip) (respectively \ip)) is represented 
with blue (respectively red) bonds. Note that Pi is applied on \ip). 
Figures nearby the diagrams are the ratios (ip\Pi\ijj) / (ip\il>) (bold) 



and e 



fo>|S i .S J -|V>)/(v#) (italic). 



b. i and j are on two distinct loops in the overlap graph of 

\tp) and \ip) and efj = 0, 

c. i and j are two sites lying at an even distance on the 
same loop in the overlap graph of \(f) and and 



'id 



-3/4. 



In the case where J,-j is uniform for nearest neighbor pairs 
of sites i and j, we introduce the nearest neighbor cou- 
pling Ji = Juj) ■ Denoting the number and net length of all 
non-trivial loops by and respectively, the number 
of trivial (i.e. length-2) loops is given by 



ni (ip,lp) -Af v ^ = 



N-L 



(10) 



For convenience we will shift away all contributions from triv- 
ial loops and rescale the Hamiltonian by a factor 4/3, and thus 
from now on we replace Hh by (4/3)"H# + J\N/2. 



In this new convention we have Eij € { — 1, 0, +1} and 



H 



<p,i> 



(11) 



where splitting the sum in Eq. (9) into two sums over trivial 
and non-trivial loops respectively yields 



Jv 



C 



<p,4> 



(i 5 j)e non- 
trivial loops 



(12) 



Because the overlap diagram of two identical configurations 



contains only trivial loops, Eq. (12) directly implies h v 



0. 



Effective Hamiltonian. Having introduced the NNVB basis 
{|V')}> the overlap matrix and the matrix H v ^, we can 
now diagonalize the Heisenberg Hamiltonian Hh projected 
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onto the singlet subspace spanned by the NNVB states. This 
could in principle be done in the NNVB basis directly, i.e. di- 
agonalizing the matrix (Wh)^.^ — {*£*YHh\4>) ■ But as we do 
not want to calculate the dual basis {|^>*)}, this might not be 
the most practicable way. Traditionally 6 one chooses the or- 
thonormal basis {\^)} = {d" 1 / 2 ^)} = {O 1 / 2 ^*)} and 
diagonalizes the matrix 

(n H ) v ^± = (<^|0 1/2 HffO- 1/2 |V) = [Hai)^, (13) 
where 

n eS = 6 1/2 h h 6- 1/2 = 6- 1/2 n6- 1/2 . (14) 

Note that on the one hand-side H e ff is expressed in terms of 
the matrices and W^^, which are both symmetric and 

easily calculated in terms of the loop structure formed by the 
NNVB states. On the other hand-side H e g arises from the 
Heisenberg Hamiltonian Hjj by a similarity transformation 
and therefore both operators are equivalent. However, work- 
ing with H e f t rather than with Hh allows for using the non- 
orthogonal NNVB basis without any explicit knowledge of its 
dual basis, which is the aim of this transformation. 
In the original scheme 617 the dual basis was omitted, but the 
definition of the effective Hamiltonian is the same in terms 
of O and %. In order to achieve this, one had to introduce 
some generalized eigenvalue problem, however missing the 
fact that H and Tin are indeed distinct operators (i.e. W is 
not the Heisenberg Hamiltonian). Furthermore, one had to 
interpret the missing * in ((p*\ip) = as emergence of so- 
called hardcore (or quantum) dimers 6 ' 7 ' 1718 . In contrast, the 
dual basis offers a quite natural framework to clarify the re- 
lations between the SU(2) (non-orthogonal) VB basis and the 
Quantum Dimer (orthogonal) basis. 



2. Operators expansion and fusions 

Diagrams. In order to explicitly calculate the effective 
Hamiltonian H e ff, we need to write down all terms in O 
and W. Noticing that the identity operator is given by 2" = 
T^x lx)(x*l> we can express the overlap operator as 

6 =ox 

= X>}<¥#)<V*I d5a) 

= ^a £ ^" 2A ^ | |, (15b) 

where Eqs. (1) and (10) have been used. As the amplitude 
V:X p (as well as /i v .^) only depends on the non-trivial loops, 
it is rather instructive to replace the flipping process \(p)(ip*\ 
by a diagram that visually describes the underlying resonat- 
ing loop structure. For example, considering the two possible 
dimer coverings around a hexagon, the resulting process is 



J 9\ 



/XX/ 



YY 



YY 



X 



(16a) 
(16b) 



where all other loops are trivial and hence omitted in the di- 
agrammatic representation. Note that this type of diagrams, 
represented in the r.h.s of Eq. (16a), actually stands for the 
sum of all flipping processes with the same amplitude and the 
same geometrical shape, i.e. all rigid motions of the diagram 
on the lattice are represented by one diagram. Furthermore, 
diagrams are symmetric, that is both processes |i y 9)('0*| and 
(ip*\ are in the same diagram. In other words we have 



(w p 0)t = u g O, 



(17) 



which we will use as determining condition for a process to 
be representable in a diagram. 

Expansion. As O and W are closely related by Eq. (11), 
a simple inspection of Eq. (1) suggests to derive the effective 
Hamiltonian defined in Eq. (14) as an expansion order by or- 
der in a. Whereas usually 617 one was to work out the overlap 
matrix in terms of the net length Cu>^> of the overlap graphs, 
we will here consider the full exponent, i.e. additionally con- 
sider the number of loops. Therefore we will denote the 
order of a graph by 2n(g) with 



n(g) = C v ^/2 - Af v ^, 



(18) 



which is clearly a non-negative integer. As an example, 

the order of the process described by Eq. (16) is 2n(gi) = 
2(6/2-1) =4. 

Then we rewrite (15) as, 



a 2 ^Co a 



(19) 



Table I compares different possible definitions of orders and 
clearly shows that the set of terms selected by truncating the 
expansion is strongly dependent of this choice. However, our 
definition might seem more natural as it actually corresponds 
to the hierarchy of occurring amplitudes in O. More impor- 
tantly, as shown in section II B 3, this choice guaranties (i) the 
locality of (i.e. non-connected processes, such as 

Aro / 

appearing in H actually disappear in T^ff at all orders) and 
(ii) the possibility of resumming all order contributions for a 
given graph. 

Using Eqs. (3) and (9) one can express the Hamiltonian as 

U = Y J W){<f\'rlH\^)m (20a) 

= Y J K^°^\v>){^\ ( 20b ) 

or, similarly as Eq. (19), 

H = ^a 2n(9) /i g w s . (21) 

9 

where h g = h v ^, for the corresponding overlap graph. 
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Processes 


Order 

2n( S ) 


Loop length 

& 


Hexagon 
number 


Degree 

dg(fl) 


/ -/ — v / / 

l$J/ 


4 


6 


1 


1 




6 


8 


1 


1 


/<> <>/ 


8 


12 


2 


2 




8 


10 


1 


1 




10 


14 


2 


2 




10 


12 


1 


1 




10 


12 


2 


1 



TABLE I: Order 2n(g) of a graph, compared to the loop length 6 and 
the number of enclosed hexagons 17 , which were defined as order in 
previous works. Note that the different expansion schemes lead to 
the very same processes but these do not appear at the same orders. 
Nevertheless, in our procedure, the leading order for a given loop of 
size C = 2p (p is an integer) in the T-i c s expansion is C — 2 (i.e. it 
scales linearly with the loop length). The last column illustrates the 
degrees of these graphs as defined in Appendix V B. 

The expression requires evaluating a non- 

integer power of O, which can be done using 

^g roTT^wW 6 - 1 '*- ,22) 

However, this still requires calculating integer powers of O 
and thus it is necessary to look at products of diagrams. Us- 
ing Eq. (17) for &\0 and &2O, it is easy to check that 
(uii^Oy = &2&iO- Therefore a simple product uiifa of 
diagrams does not fulfill Eq. (17) and hence cannot be repre- 
sented as a diagram. On the contrary, the symmetrized form 

u 9 = - {wi,a> 2 } (23) 

obviously verifies Eq. (17) and can be represented as a graph. 

Fusions. The issue of generating graphs out of symmetrized 
products of graphs is called fusion of diagrams. These fusion 
rules can be derived straightforwardly, by using the definition 
of the participating processes, such as Eq. (16). In a first step, 
one will only derive fusions of terms appearing in the expan- 
sions of the Hamiltonian (Eq. 21) and in the overlap matrix 
(Eq. 19). These processes are always flipping a dimer config- 
uration on the boundary of the representing diagram and are 
therefore called kinetic (or off-diagonal) terms. But generally, 



a fusion can create processes which do not appear in the over- 
lap matrix nor in the Hamiltonian. However, Eq. (23) guaran- 
tees, that there is a diagrammatic representation for every of 
such fused diagrams. 

Let us illustrate in details this issue by discussing the fol- 
lowing fusion rules, that are valid on the kagome lattice. 



^{/n/,/n/} =/n/ + 1 /rrrr / (24a) 




The first example, see Eq. (24a), shows that the fusion of 
two identical kinetic diagrams always involves the emergence 
of potential (or diagonal) terms that are drawn here with a 
yellow shape. This kind of terms arises when a flipping pro- 
cess is acting twice at exactly the same position. Although a 
repeated flip leaves the dimer configuration unchanged on the 
plaquette, it is important to emphasize that it is not identical to 
the identity operator as it is actually checked whether the pla- 
quette is flippable or not. Indeed, applying a flipping process 
to a non-flippable plaquette annihilates the state. This can also 
be seen from the fact that contributions of potential terms can 
be expressed as \ip)(ip*\ and as an example on a hexagon we 
have 



^92 = /<T 7 (25a) 

= i /^7)( /^7*i + 1 /o7)( /^7*\- < 25b > 

Furthermore, the fusion of diagrams always generates dis- 
connected diagrams, where the flipping or checking processes 
are simply happening at distant positions (e.g. last terms on 
the r.h.s of Eqs. (24)). Notice, that due to combinatorics the 
prefactor is not the same in the case when identical or different 
diagrams are fusing. 

Equations (24b) and (24c) illustrate, that when both kinetic 
or potential processes are happening close to each other (i.e. 
sharing at least one bond), the resulting diagrams generally 
look more complicated. A typical merging of two loops into 
one larger loop is displayed by Eq. (24b). Furthermore, more 
unusual assisted kinetic processes where the plaquette flipping 
requires the presence of specific trivial loops in the neighbor- 
hood of the plaquette, see Eq. (24c), also emerge from the 
fusion rules. Notice that the kinetic and potential processes 
are identical on trivial loops and correspond to verifying the 
presence of dimers at a specific position. 

While the details of the fusions rules are lattice-specific, 
the key properties of fusions for the derivation of an effective 
Hamiltonian are actually quite general. In Appendix V E, we 
introduce a general lattice-independent diagrammatic notation 
in which each connected part of a (disconnected) diagram 
((#) (#) . . . (#)) is generically represented as (#). General 
fusion rules (see Appendix V F) produce new connected terms 
by connecting two or more connected parts. Such terms are 
generically represented as (•-•) (fusion of two connected 
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parts). For example, all the rules given in Eqs. (24) lie in the 
generic class of rules, 



((•) x (•)>=<(•"•)> + ((• •)). 



(26) 



Note that the notation (.) defined in Appendix VE, conve- 
niently absorbs all combinatorial prefactors. Fusing more 

connected parts leads to diagrams such as (((-^~))) or 
((O-O-O)) (fusions of three connected parts). Let us remark 
that while both terms are connected, these two terms are in- 
equivalent since in the latter case the first and third part do not 
fuse. 

Provided these general fusion rules, one can work out the 
effective Hamiltonian H e ff in a systematic way. It is important 
to note that O and ft do not contain any potential term. They 
only emerge by fusing kinetic processes. On the other hand, O 
and ft contain disconnected kinetic terms and when fusing di- 
agrams we always generate those non-local processes as well. 
Such a non-locality would be physically inconsistent in the re- 
sulting effective Hamiltonian fL eS = 0~ 1/2 fLO~ 1/2 . How- 
ever, all the non-local terms disappear while deriving ft e ff and 
thus the effective Hamiltonian is local as the initial Heisen- 
berg Hamiltonian. As this property, proofed in the next sec- 
tion, is valid at all orders in a, this holds for the exact all-order 
resumed effective Hamiltonian as well as for any truncated 

As- 



3. Proof of locality and re summation scheme 

Locality. Comparing equations Eqs. (19) and (21), the only 
difference between O and ft is the additional weight h g for 
the Hamiltonian. Therefore it is convenient to define a gen- 
erating function Z(a, (3, /i) (defined in Eq. (57) in Appendix 
V B), such that both O and ft can be obtained from it by 



0=1 



O = Z]^, 
0=1 



(27a) 
(27b) 



where /? and /i are some internal parameters (for details see 
Appendix V B). Therefore, ft e ff can be recast as 



0=1 



(28) 



After some standard operator manipulations (see Appendix 
V C), one can show, that the last equation is equivalent to 



H 



eft 



\ " 



0=1 P=0 



1 ( 


r - \ 2 p 
InZ,) 


2 2p ( 


2p+l)\ 



dAnZ. 



(29) 



involving only iterated commutators of In Z on In Z . 

Now, supposing that In Z only contains connected terms, 
it is clear that this cannot be different in c^ln Z . Moreover, 



the commutator of two connected diagrams cannot be discon- 
nected: (i) either the two diagrams act on disconnected pla- 
quettes, and hence commute, (ii) or they lie on neighboring 
plaquettes (i.e. sharing at least one bond) and, according to 
the concept of fusion introduced in the previous section, they 
fuse to a single connected diagram. This can also be shown 
easily using [A, [A, B]} = {B, {A, A}} - {A, {A, B}} and 
the diagrammatic scheme in Appendix V E. 

It turns out, that In Z indeed only contains connected terms, 
which is an important result both from a conceptual and aes- 
thetic points of view: it shows that "Heft is local and therefore 
provides firm grounds to the consistency of the method. Con- 
versely, producing an effective Hamiltonian with non-local 
terms would put the whole scheme into question. 

The quite technical demonstration of this point is split in 
several Appendices V D to V G. Let us sketch the various steps 
of the proof : 

1. We establish a linked cluster theorem, i.e. express the 
logarithm of the generating function as sum of cumu- 
lants (see Appendix V D). 

2. We set up a lattice independent diagrammatic notation 
(Appendix V E) and work out general fusion rules (Ap- 
pendix V F). 

3. Using these rules, we establish Eq. (99) showing that 
every cumulant can be reexpressed as a combination of 
lower order cumulants. This allows finally to show by 
mathematical induction that cumulants represent con- 
nected processes and therefore \nZ does not contain 
non-connected terms (see Appendix VG). 

As we will see, this generating function formalism is not 
only useful to establish the locality of the effective Hamil- 
tonian but provides also a practical framework to resum the 
series in a, giving the amplitudes of dominant terms in ft e e. 

Fully connected diagrams. A special class of connected di- 
agrams, called fully connected diagrams is produced by com- 
plete fusion of all parts (i.e every diagram is connected to ev- 
ery other diagram by sharing at least one bond). In the lattice- 
independent diagrammatic representation, we denote such di- 
agrams as 




These diagrams are very important for two reasons: 

• This class of diagrams includes the most compact terms 
(i.e. involving short loops) of the effective Hamilto- 
nian that are obtained by fusing identically shaped di- 
agrams exactly at the same place (e.g. the first term 
in the r.h.s of Eq. (24a)). These processes are also the 
dominant terms in the a-expansion of "H e ff- Indeed, 
the leading order for a length-£ connected process in 
<U eS = G-^UG- 1 / 2 is a c - 2 . 
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Their exact weight in In Z is particularly simple to ob- 
tain, as explained in Appendix V H and, 



(-iy 



rn 



m— 1 

less connected terms. 




(30) 



Elementary diagrams & Resummation. We define the ele- 
mentary diagrams as the set of processes (kinetic or potential) 
that cannot be cut into two sub-processes. For example, 



/07 and 



are elementary while 




(31) 



(32) 



is not since it can be cut into the two parts of Eq. (31). Obvi- 
ously, all elementary graphs appear as fully connected graphs 
in In Z: for example uj g2 given by Eq. (25) is a full connection 
of u) gx , see Eq. (16), with itself. Notice that there exist fully 
connected terms that are not elementary (e.g. Eq. (32) is a full 
connection of Eq. (31), but is not elementary). 

Generic elementary kinetic and potential diagrams will be 
represented respectively as, 



© and © 



(33) 



Fusion rules of elementary diagrams, of which Eq. (24a) is a 
typical example, can be compactly written as 



\{m®}=w 



0© 



©© 



©© 



(34a) 
(34b) 
(34c) 



In order to compute the full weight of elementary diagrams 
in In Z, we need to track all occurrences of these terms in 
Eq. (30). High order contractions of these processes are typi- 
cally given by iterating Eq. (34). But obviously, the last (dis- 
connected) terms of Eq. (34) cannot produce fully connected 
terms by further iterations. As a consequence, the evaluation 
of In Z only requires the simple reduced rule, 



© for odd n, 



(35) 



©] for 



even n. 



Since a length-£ elementary process with weight h in the 
expansion of H occurs at order a c ~ 2 , its contribution to the 
generating function is, 



/3a 



C-2 e vh 



© 



(36) 



Using Eq. (30), the relevant contribution to In Z is 
In Z = . 



(_i ) m 

E/3 a m(C-2) c m^h 
m 



which is easily evaluated, using Eq. (35), as 
1, 1 + I3a c - 2 e^ h 



\nZ = 



In 
2 1 



© 



(37) 



(38) 



+ iln(l-/3 2 a 2 ( £ - 2 )e 2 ^) [©]+... 

Evaluating H e ff is done using Eq. (29) which requires in 
turn to compute c^ln Z and all its iterated commutators with 
InZ. But since, by definition, elementary diagrams cannot be 
produced by fusing smaller diagrams, the only relevant terms 
to the amplitude of Eq. (33) in H e ff produced by orders p > 
in Eq. (29) are (i) commutators of an elementary diagram with 
itself at the same place or (ii) commutators of two elementary 
diagrams acting on disconnected places. Obviously such con- 
tributions vanish identically. As a direct consequence, only 
order p = of Eq. (29) contribute to the amplitude of ele- 
mentary diagrams in H e g. In other words, for any elementary 
process 



■Heff = dJnZ. 

0=1 



(39) 



It is then straightforward, using Eq. (38), to evaluate the con- 
tribution of elementary diagrams to H e ff : 



"Heff — 



o 



£-2 



1 - a 2 ( £ - 2 ) 
q 2 ( £ - 2 ) 

1 _ a 2(£-2) 



(40) 



© 



where C is the length (see Table I) of the elementary (kinetic) 
diagram and h its "bare" energy as appearing in the expansion 
of H, given in Eq. (21). 



III. GENERALIZED QUANTUM DIMER MODEL FOR 
THE KAGOME ANTIFERROMAGNET 

In this section we will apply our scheme to the kagome lat- 
tice and consider the nearest neighbor coupling Jj = J > 
only. First, we make a brief review of the current understand- 
ing of this model. 



A. Current status 

The properties of the QHAF on the Kagome lattice have 
been actively explored for the last two decades by analytical 
and numerical techniques. Although a number of important 
results have been obtained, the nature of the groundstate is 
still a mystery. Lanczos ED of small clusters 19,20 have clearly 
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Processes 





O 



U/J 



-3a 4 



-2a b 



-2a b 



-2a b 



-6a a 



Processes 




O H/J 



-5a 1 



-4a 1 



a 10 -4a 10 



TABLE II: Expansion of X — O — 1 (or 6 excluding order 0) and 
H up to order 2n(g) — 10. 



Processes 



/Q7 



4£X. 



/£// 



4t£ 



Weff/ J 



LO 



-3a 4 



3a s 



-2a b 



-2a b 



-2a b 



2a 1 



2a 1 



2a 1 



3n* 



2a° 



2a 



Processes 




'Hat/ J 



LO oo 



n- 8 « 8 .,. 



1-a 16 



J__ a 16 



1-a 16 



revealed an exponentially large number (w.r.t. system size) 
of singlets below the first triplet excitation (spin gap), a fact 
reminiscent of the NNVB basis 12 ' 21 . The accessible cluster 
sizes remain however too small to definitely conclude whether 
the spin gap survives in the thermodynamic limit 22 . Recent 
DMRG studies 23 seem however to support it. 

A number of studies have also been devoted more pre- 
cisely to the nature of the groundstate itself. For exam- 
ple, an early ED analysis of the four-spin correlations (i.e. 
dimer-dimer) first pointed towards a short-range dimer liquid 
phase 24 . Alternatively, translation-symmetry breaking VBC, 
two-dimensional analogs of the Majumdar-Gosh 25 dimerized 
chain, have been proposed on the basis of various other ap- 
proaches like large-N SU(N) techniques 4 ' 26 and mappings 
to low-energy effective hamiltonians within the singlet sub- 
space 27-29 . Recent series expansions around the dimer limit 30 
showed that a 36-site unit cell (i.e. a 2\/3 x 2\/i super- 
cell of up-triangles) VBC would be preferred in agreement 
with Refs. 26,29. However, the interpretation of the ED low- 
energy singlet spectrum 31 on the 36-site periodic cluster re- 
mains problematic not showing the expected quasi-degenerate 
VBC groundstates. Therefore, it has been proposed 22 that the 
Heisenberg model might be within (or in the close vicinity of) 
a spin liquid phase such as the algebraic spin liquid 5 . We now 
turn to the application of our procedure which offers a new 
versatile scheme to investigate non-magnetic singlet phases 
as VBC. 



TABLE III: Elementary processes in H e ff at Leading Order (LO - 
see text for explanation) and fully resummed series (oo - see Eq. 40). 



B. Expansion 



Leading orders. Table II shows all the contributing terms 
in O and % up to order 2n(g) = 10. Of course, as VtV = 1 
and T-L VtV = 0, the identity process 1, which only contains 
trivial loops, appears with order 2n(g) = in O but does not 
contribute to H. Therefore we define O = 1 + X, where 
A" is a short hand notation for all processes in O, except the 
identity. 

Using equation Eq. (22) the effective Hamiltonian "H ett - = 
O- 1 / 2 ^- 1 / 2 can be written as 



H 



eff 



- % -2 



\ {*,%} (41) 

\{*,{*,*}}+\{M*>*}} 

terms with at least 3 anticommutators. 



Thus we can immediately see, that all elementary kinetic pro- 
cesses have, at the lowest order, the same weight in % and in 
Heff (compare Table II with III). Additionally considering the 
second term in Eq. (41), one can confirm the amplitude of 



to be — (a 4 ) ■ (— 3 Ja 4 ) = 3 Ja 8 . Here, the first term is due 
to the fact that H does not contain any potential graph (see 
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Table III), whereas the second term illustrates the fusion of 
kinetic diagrams in X and H to potential ones (see Eq. 24a). 
Similarly, it is easy to see that the amplitude of 




cancels out to (-6Ja 8 ) - 2 • (a 4 ) ■ (-3Ja 4 ) = 0, which 
is a direct verification of the general property proved in Sec- 
tion II B 3 stating that disconnected graphs do not contribute 
to the effective Hamiltonian at any order. The leading orders 
of all the diagrams presented in Table III can be obtained con- 
veniently using a direct evaluation based on Eq. (41). 

Higher orders. Diagrams with larger leading orders 
(typically enclosing several hexagons) as well has further 
renormalizing corrections to low leading order diagrams 
are more problematic to obtain using simple arguments. 
Nevertheless, the calculation can be systematically extended 
to significantly higher orders using the expansion scheme 
presented in Section II B 2. Up to the chosen order, this 
requires (i) the enumeration of all terms appearing in T-L, (ii) 
a careful enumeration of fusion rules that proliferate as the 
order increases, (iii) the expansion of 0~ 1/2 and (iv) the 



evaluation of O^^HO^ 1 / 2 . Note that the two last steps 
explicitly require using the fusion rules obtained in step (ii). 
The details and results of such a procedure up to order a 14 
are too lengthy to be presented in this article and are therefore 
provided as supplementary material at [URL will be inserted 
by AIP] in which we include all the relevant fusion rules up 
to order a 14 as well as extensions of Tables II and III. 



C. Resummation 

A simple inspection of Table III reveals that all the lead- 
ing processes in the effective Hamiltonian are elementary di- 
agrams in the sense defined in Section II B 3. Indeed, none 
of the terms enclosing only one hexagon can be split in sub- 
processes. This important remark shows that the resummation 
scheme of elementary diagrams presented above applies di- 
rectly. In particular the equation (40) immediately leads to the 
resummed amplitudes of the effective Hamiltonian presented 
in Table III. The explicit form of % e ff is obtained by setting 
a = ij \/2 which leads to 



D. Discussion 

General remarks. Note that the kinetic part of this Hamilto- 
nian 32 ' 33 is quite close to the one originally proposed by Zeng 
& Elser 17 with only small differences in the magnitudes of 
the processes, differences introduced by our infinite order re- 
summation scheme. This provides strong evidence that the 
expansion indeed converges rapidly. However, very impor- 
tantly, our Hamiltonian includes also diagonal (i.e potential) 
terms which turn out to play a major role but which were 
not included in equation (7) of Ref. 17 based on a different 
expansion scheme 34 . Indeed, the low-energy gap presented 
on figure 3 in Ref. 17, whose magnitude is approximately 
36 x 0.0055J ~ J/5, splits two sectors with 2 and 1 flip- 
pable hexagon(s) respectively (see figure 4 in the same ref- 
erence). Interestingly, the value J/5 is precisely the ampli- 
tude of the potential term disadvantaging flippable hexagons 
present in the effective Hamiltonian (42). This strongly sug- 
gests that this gap, known to be absent from the exact spec- 



I 

trum of the Heisenberg Hamiltonian 35 , is an artifact produced 
by truncating the expansion. As already mentioned (see Sec- 
tion II B 2 and Table I) the traditional expansion scheme 617 , 
by modifying the hierarchy of terms in the expansion, as a 
tendency to push away to higher orders in a the emergence of 
processes in the effective Hamiltonian. On the contrary, in- 
cluding such a potential term is likely to close this gap and 
brings the Hamiltonian given in Eq. (42) closer to the actual 
gapless low-energy phenomenology of the kagome antiferro- 
magnet. 

Another important remark is that the amplitude of the ki- 
netic (and potential) pinwheel process, 




denoted Jn below, exactly vanishes at all orders. As dis- 
cussed below including a finite J12 in the model lifts a very 
special degeneracy of the groundstate manifold. 
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Large scale numerical computation. In contrast to the case 
of the frustrated square lattice 36 , for the kagome lattice we 
obtain an effective model whose leading coefficients have al- 
ternating signs precluding any stochastic approach. However, 
Lanczos exact diagonalizations can be preformed on relatively 
large clusters 7 due to the very constrained nature of the dimer 
basis that greatly limits the number of states (2 3 + 1 compared 
to 2 N for SU(2) spin-1/2 models). Furthermore, group the- 
ory techniques can be applied to block-diagonalize the Hamil- 
tonian matrix in each of its irreducible representations (IR- 
REP) hence further reducing the practical number of degrees 
of freedom. For example, for the most interesting clusters 
with N = in 2 or 9m 2 sites (which possess all relevant space 
group symmetries of the infinite lattice) such as the 36-, 48-, 
108- and 144-sites clusters, the increasing Hilbert space sizes 
of their smallest (largest) IRREP are roughly of the order of 15 
(170), 70 (2xl0 3 ), 80xl0 6 (950xl0 6 ) , 200xl0 9 (3xl0 12 ) 
respectively. Hence, current supercomputers enable to tackle 
the 108-sites cluster while the larger 144-site cluster might be 
reachable within a few decades. 

Model properties. The numerical results for the Hamil- 
tonian given in Eq. (42) presented in Ref. 7 (corresponding 
in fact to an earlier extremely close 14-th order estimation 
of the model) have been summarized in the phase diagram 
of Fig. (3). Here, we have introduced two extra parameters; 
(i) A to interpolate (linearly) between the "Heisenberg point" 
(see Eq. (42)) at A = 1 and the "RK-point" of Ref. 37 (de- 
noted T-Lrk below) at A = 0. (ii) A finite pinwheel amplitude 
J12. This double interpolation can be summarized in the two- 
parameter Hamiltonian 

Hmterp. (7, Jvi ) = T^eff + ( 1 - l)H R K (43) 



■A2 + 4(7-1) 



where J = 1 and the RK-Hamiltonian has been defined in 
Ref. 37 withT = -1/4 32 . 

In agreement with series expansions 30 and earlier work 26 ' 29 
on the QHAF, the GS of our corresponding effective model 
(see Eq. (42)) is found to be a VBC with a 2^ x 2\/3 su- 
percell of up-triangles (carrying 36 sites) whose underlying 
dimerization pattern corresponds to an honeycomb-lattice ar- 
rangement of the resonating "perfect" hexagons (colored in 
blue in Fig. (3)). Because J12 = 0, the Heisenberg effective 
model given in Eq.(42) is also characterized by a degeneracy 
of the even and odd resonating pinwheels of the motive (col- 
ored in yellow in Fig. (3)). As pointed out in Ref. 30 this leads 
to an extra Ising-like degeneracy. Remarkably, the numerical 
results show also that the "Heisenberg point" lies very close 
to the critical line which separate the VBC phase from an ex- 
tended dimer liquid phase similar to the one at the "RK point". 
Interestingly, a general field-theoretic framework 38 based on a 
double Chern-Simons theory correctly describes such a Quan- 
tum Critical Point: one considers the spectrum of visons (i.e. 
topological defects) in the dimer (so-called Z2) liquid phase, 
and studies how they condense - the condensation of visons 
leads to VBC order. Our approach applied to the Kagome lat- 



VBC 36-site 





QCP I 1 **X 
I A 
l 

VBC 36-site 
"odd" 




FIG. 3: (Color online). Semi-quantitative phase diagram of the ex- 
tended GQDM for the Kagome lattice as a function of the pinwheel 
resonance amplitude J12 and the parameter 7 (see Eq. (43)). The two 
dashed blue lines represent (i) an interpolation AH e ff + (1 — X)Hrk 
between the RK model and the effective model for the Heisenberg 
quantum antiferromagnet and (ii) a cut at 7 = 1 along J12, and cor- 
respond to the simulations performed in reference 7. Red lines are 
qualitative phase transitions. Note that for 7 > 0.9 a finite J\i lifts 
the degeneracy between two VBC's with identical 36-site unit cells 
but opposite parities P± of their resonating pinwheels (in yellow). 



tice as well as other numerical studies of a generic QDM on 
the triangular lattice 39 strongly support such a scenario. Even 
though the Z2 phase has no broken symmetry, it does have a 
topological order. In principle, the topological order can co- 
exist with the VBC, and so their disappearance at a common 
critical point without fine-tuning, can be considered as a non- 
LGW (Landau-Ginzburg-Wilson) transition 40 . 



IV. CONCLUDING REMARKS 

In summary, a systematic non-perturbative method has been 
developed to describe quantitatively the low-energy physics of 
frustrated QHAF. Provided that the latter is governed by fluc- 
tuations of (short-range) singlets (which is believed to happen 
in many cases), this method is fairly general and can be ap- 
plied in principle to any lattice geometry. The low-energy ef- 



11 



fective model takes the form of a Generalized Quantum Dimer 
Model which is proven to be local, a physical requirement. 
Complete practical formalism is presented as well as other 
important general results. It is also shown that the expansion 
scheme can be pursued up to high orders and the most relevant 
terms can be resummed. 

As a practical implementation of the method, we have con- 
sidered the much debated Kagome QHAF where a fully re- 
summed parameter-free GQDM can be obtained. Although, 
the resulting model bears a sign problem (e.g. could not be 
simulated by QMC techniques), ED results up to clusters with 
108 sites can be performed 7 showing strong evidence in fa- 
vor of a large supercell VBC. Also, it turns out that, in some 
extended parameter space, the effective model of the Kagome 
antiferromagnet lies in the close vicinity of a critical line to- 
wards a topological quantum dimer liquid phase, somehow 
clarifying the low-energy puzzle of the original spin model. 
Let us recall that a similar approach has also been applied 
to the frustrated square lattice 36 . This shows that effective 
GQDM can efficiently describe the most frustrated quantum 
magnets and greatly help to understand their properties on 
larger (smaller) length (energy) scales compared e.g. to stan- 
dard ED techniques. It could be used also in 3 dimensions to 
tackle the quantum antiferromagnet e.g. on the hyper-kagome 
lattice. 

A number of Hamiltonian extensions could easily be in- 
cluded in the following approach like (i) other SU(2)-invariant 
terms as multiple exchange or longer-range exchange interac- 
tions 44 , (ii) doping with static or mobile holes 42 and (iii) the 
inclusion of spinons 43 and triplets. In the case of the Kagome 
lattice, it is important whether small non-SU(2) terms like 
Dzyaloshinski-Moriya interaction could also be inserted in a 
realistic way. All these issues are left for future studies. 

This work is supported by the French ANR program ANR- 
08-JCJC-0056-01. 



V. APPENDICES 
A. Dimers, overlaps and sign conventions 

In this appendix we recall some of the most important basic 
properties of VB states and two orientation conventions that 
are suitable for the derivation scheme proposed in this article, 
respectively for bipartite and non-bipartite lattices. 

Overlaps. Let us consider a Valence Bond (VB) wavefunc- 
tion 



\<f) 



N/2 



fc=i 



where 



V2 



(ttiW - \UW) 



(44) 



(45) 



is the SU(2) dimer wavefunction on sites i and j. Clearly, this 
dimer is antisymmetric with respect to the permutation of i 
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\ 
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FIG. 4: Explicit bosonic (e = +1) and fermionic (e = — 1) conven- 
tions for the 15 VB states on 6 sites. 



and j, which makes it necessary to use arrows in the graphi- 
cal representation (see Fig. 1), in order to distinguish between 
all the 2 N I 2 possible orientations of a dimer covering. In the 
following lines, we emphasize two particularly useful orien- 
tation conventions. This allows to omit arrows on the dimers, 
provided we give a recipe to keep track of the relative signs of 
the VB states given by Eq. (44). Two arbitrary VB states \ip) 
and \tp) are non-orthogonal and 



(46) 



where ni is the number of loops in the overlap diagram (see 
Fig. 1) and r) (ip, i[>) is a sign that depends on the chosen con- 
vention for orientating SU(2) dimers. 

The question of orientating dimers in order to fix the form 
of r) v ,^ is determined in turn by two parameters: (i) the nature 
(bipartite or non-bipartite) of the lattice on which the dimers 
are constructed and (ii) the constraints one puts on the type 
of dimers considered (arbitrary long range dimers, arbitrary 
long range dimers with a bipartite constraint, nearest neighbor 
dimers,. . . ). In the perspective of performing an expansion of 
the overlap matrix as powers of a small parameter it is essen- 
tial to choose conventions where the sign in the power law of 
equation (46) can be systematically absorbed: 



(¥#) 



v iV-2n i ( ¥ >,V) 



(47) 



The only two possible conventions denoted b (bosonic, +) and 
/ (fermionic, — ) for which 



(48) 
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correspond to «j = 1/ \/2 and at = i/ \[2. These two rep- 
resentations, widely discussed in the literature (see for exam- 
ple Ref. 41), originate from the fact that dimer wavefunctions 
can be constructed out of a vacuum by applying bosonic or 
fermionic operators resulting in different statistic and hence 
sign conventions. In the following paragraph we will fully 
choose to present these two conventions without any direct 
reference to bosonic or fermionic operators, emphasizing the 
fact it can be directly understood as a conventional orientation 
of dimers, as made explicit in Figure 4 for 6 sites. 

Bosonic convention on a bipartite lattice. In this conven- 
tion, the set of N sites is divided into two iV/2-site subsets 
A and B, and dimers are oriented such as i 6 A and 
j € B. As a direct consequence of this orientation, arrows are 
diverging from all A sites and converging to all B sites in the 
overlap diagram. The contribution of each loop is thus +2 and 
a = 1/V2. 

Note that this convention can also be used on non-bipartite 
lattices, provided a (arbitrary) bond orientation is specified. 
However, this generally produces a sign in the overlaps that 
cannot be absorbed such as in Eq. (47). 

Fermionic convention. Choosing an arbitrary numbering of 
sites, we start from a reference configuration, 



|^o> = [1, 2] [3, 4] ... [AT — 1, N}. 



(49) 



Here on a non-bipartite lattice, in the bosonic representation of 
the dimers we need to specify the sign of each dimer following 
some prescription. Therefore, by convention, we write the 
dimer wave function defined in Eq. (45) such as i/. < jk- In 
the expression 



N/2 

n [ikjki 
fc=i 



(50) 



the symbol * means that < jf.. Any state (50) can be ob- 
tained from \(po) by applying a unique permutation it acting 
on site numbers which generically writes, 




N-l N 

lN/2 j/AT/2 



(51) 



Then we define 



N/2 

\<P) = Cr(71>) Y[ [ik,3k] 
fc=l 



0-^)71-^0), 



(52) 



where cr(ir) is the signature of the permutation ir : 

i<j 



a{-K) = n 



3 



(53) 



Note that the state defined in Eq. (52)) is in fact independent 
of the prescription to order the individual dimers but depends 
only on the ordering of the sites in the reference state given 
by Eq. (49). An explicit example of this convention (as well 
as the bosonic convention) is given on figure 1 . Let us remark 



that if \(p) and \ip) denote two states deduced from the refer- 
ence state given in Eq. (49) by applying the permutations ir v 
and ir^f, then their relative sign according to Eq. (52) is <j(tt) 
where ir = n v n^,. 

We are going to show that the convention of Eq. (52) in- 
deed implies r){tp,ip) = (-i)Mv,Tp)-N/2 in Eq _ (46)) _ L et 

us compute (ip\Pi\x) where Pi is the permutation operator on 
the bond I = (i, j). As presented on figure 2, three cases can 
occur: 

a. The bond I connects two sites lying at an odd distance 
on the same loop. The total number of loops of the 
overlap diagram remains unchanged and a sign change 
occurs: (tp\Pi\x) = -{<p\x)- 

b. The bond I connects two distinct loops. Two loops are 
merged into a single one in the overlap diagram and no 
sign change occurs: (ip\Pi\x) = +l/2((p\x). 

c. The bond I connects two sites lying at an even distance 
on the same loop. One loop is created in the overlap di- 
agram by slicing one loop into two and a no sign change 
occurs: (f\Pi\x) = +%{<p\x)- 

Obviously, the state Pi\x) violates the convention of 
Eq. (52) since it is deduced by applying a 2-site permutation, 
denoted 7r;, on a conventional state without being multiplied 
by the factor (7(717) = —1 (note that ).The corresponding con- 
ventional state is thus — o~(tti)Pi\x) an d the three cases 
considered above become: 

«• (¥#) = (¥>|X>- 

b. (<p\i>) = (-1/2)MX>. 

c. = (-2)Mx). 

The group structure of the permutations ensures that (i) the 
signs of all states defined by Eq. (52) can be consistently gen- 
erated from the reference state given by Eq. (49) and (ii) a 
change An; in the number of loops comes with a sign change 
(-l) An '. This indeed implies 77 (tp, tp) = (-i)M<p,i>)-*/2. 



B. Operators and Generating function 

In this appendix, we introduce the generating function Z. 
As explained in Section II B 2, the Hamiltonian and Overlap 
operators can be expanded in powers of a as linear superposi- 
tions of dimer flipping processes uj g : 



(54a) 
(54b) 



Here 2n(g) denotes the order of the graph g representing 
the process. Note that ui g is a shortcut notation for an implicit 
summation over all the lattice of all plaquettes with the shape 
g. By definition ujq = 1 and ho = 0. 
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The degree dg(<?) of a graph g is the number of its con- 
nected parts. By definition, if a) is a connected graph then 
dg(u)) = 1. By convention, dg(l) = where 1 is the identity 
operator. For some examples, see Table I. 

By introducing the X g operators defined as, 



X g = e» h °u, g , 



(55) 



such as d^Xgl^Q = h g Cj g and X g \^Q — tu g , one can ex- 
press H and O as 



0=1 

6 = z\^ =0 , 

0=1 



(56a) 
(56b) 



using the generating operator Z defined as, 

Z(a,(3^) = Z 2 aM9) P M9) X 9 - (57) 



In turn, the effective Hamiltonian, 

•Hetf = 6- 1/2 h6- 1/2 , 

is rewritten as, 

n efi = q z- 1 / 2 (d,z) z- 1 ' 2 . 
0=1 



(58) 



(59) 



C. Hat as a function of In Z 

The purpose of this Appendix is to explicit the transforma- 
tion by which given by Eq. (59) is rewritten using only 
In Z, as announced in Eq. (29) of Section II B 3. 

We use the general result, 

d(e F ) = t e^ F (dF)e sF ds, (60) 
Jo 



for F = In Z to express H e ff as, 



m=o J_i 

0=1 2 



(^lnZ) 



Z s ds. 



(61) 



We now use the Hadamard formula 



e A Be~ A = B = B + [A, B] + ~ [A, [A, B]] + ... (62) 
with A — s In Z and B = 9^ In Z to rewrite Eq. (61) as, 



"Heff = | / < 
u=0 
/3=1 V~ 5 



InZ. 



ds ajnZ. (63) 



Other equivalent forms are, 



Heff = 2- 
ju=0 
(9=1 



sinh f 1 


InZ,) 







c^lnZ 



(64) 



or, by expanding sinh into a Maclaurin series, 



fi=0 ' — ' 

0=1 p=o 



1 ( 


r - \ 2 p 

InZ, J 


2 2p ( 


2p + l)! 



a„inz. 



(65) 



D. Cumulants 

In this appendix we derive the expression of In Z enter- 
ing in Eqs. (63), (64) and (65) as a (3 expansion of non- 
commutative cumulants. 

Let us introduce the collection of all a 2n ^X g of same de- 
gree d: 



a 2n ^x g . 

g such as 

dg(g)=d 



(66) 



z^ 1 ' involves all simply connected graphs, z^ all graphs of 
degree 2, etc and £W = i. Combining Eqs. (57) and (66) 
leads to 



Z = Y,f3 n z (n \ 

n=0 



and 



In Z 



E 



1 / z -l 



(67) 



(68) 



Let us call 5™ the set of all sequences s = {sj} made 
of strictly positive integers Sj such as the number of integers 
\s\ = n and J^ies Si = m - Notice that considering non- 
vanishing integer sequences, the minimal value of J^ies s i f° r 
s £ S™ is n and therefore S™ is non-empty only for m > n. 
Moreover, S™ is the empty set for all m. Using this notation, 
one can regroup all terms with J^ies Si ~ m ' anc ^ therefore 
we obtain 



InZ = - 



(-1) 



oo oo 

E E n2 

n— 1 m—n s£<S m \i=l 



LfrY, II{^ i, (69) 



which, by inverting summations, can be written in turn as, 

oo 

Ini = £ /3 m zi m \ (70) 

m— 1 

where the order-m non-commutative cumulant zi" 1 ^ is defined 
as. 



(-1)' 



^ n2 

n—1 s Go "* \1= 1 



e n{^ U (7i) 
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for m > 0. By definition zi°^ = 1 . 



The explicit expressions of the first zi™' are: 



(72) 



( 3) =^ 3) -i{5 (1) ,i (2) } + ^{5 (1) ; {5 (1) ,i (1) }} 



2 (4) = £ (4) _ 



Interestingly, these expressions of z^ in terms of z^ can 
be inverted using the fact Z = cxp(ln(Z)). In the spirit of 
Eqs. (68) and (69), 



oo 1 



m=0 
oo oo 



(73) 



EE 

ml 2 m 

e n{ 

771—0 n—m s£S n 



n=0 m=0 \i=l / 



This last expression allows to identify, 



n / n \ 

m— 1 s€$Il / 



Here is the explicit form of the first terms, 



5 (D = 



(75) 



^ ) = 4v + l{4 1 \^} + y 4 {4 1 \{4 1 \4 1) }} 
+ ^K>,{^,jf>}} + l{if>,{e4 l >}} 



£. Diagrammatic notation 

This Appendix is devoted to the introduction of a lattice- 
independent diagrammatic notation. Let us introduce the fol- 
lowing diagrammatic notation : 



(0 0-0)=n!z< n \ 



(76) 



and its cumulant counterpart, 



(••-•) =m\zi m \ 



(77) 



Defining, for any s — {si, . . . , s m } € 5^, the combinato- 
rial factor 



i (Epli s pj ! 



(78) 



we rewrite Eqs. (71) and (74) as, 
((• •-•),) = 



(79) 



£(-l)«-i(n - 1)! £ [*]<(• •-•) x • •• x (• •-•)) 



ses:; 



((OO-O)) = (80) 

n 

n 

x; e [^(••••••) c x •••*(••••••),>, 



m— 1 s£5" 



where the brackets (.) notation stands for: 



(81) 



The diagrammatic counterparts of Eqs. (72) and (75) are 
respectively 



(82) 



((•)c> 


= 


[!]((•)) 


( 


<(o o) c ) 


= 


[2]((« •)) 






- 1! 


[l,l]((«)x(»)> 




((• • •)„) 


= 


[3]((« • •)) 






- 1! 


[l,2]((«)x(» •)) 






- 1! 


[2,1]((« »)x(»)) 






+ 2! 


[l,l,l]((«)x(«)x(«)) 




<(o o o o) c ) 


= 


[4]((0 0)) 






- 1! 


[l,3]<(«)x (•••)) 






- 1! 


[3,l]<(#»#)x(«)> 






- 1! 


[2,2]<(« m)xcm)) 






+ 2! 


[U,2]((i)x(#)x(i •)) 






+ 2! 


[l,2,l]((«)x(« «)x(«)) 






+ 2! 


[2,1,1]((« • x (•) x (•)) 






- 3! 


[l,l,l,l]((«)x(«)x(»)x 


(•)) 



15 



and 

((•)) = [1] ((•).) (83) 
((• •)> = [2]<(« •)„> 

+ [l,l]<(») c x(«) c ) 
((•••)) = [3] ((OOO)J 

+ [i,2]((«) c x(» m) c ) 

+ [2,1]<(« •) C X(«) C > 

+ [i,i,i]((») c x («) c x (m) c ) 

((0 0)) = [4]((0 0) c > 

+ [l,3]((«) c x (•••)„) 

+ [3,1] ((•••), x (m) c ) 

+ [2,2]<(« •),. x (• •),.) 

+ [i,i,2]<(«) c x («) c x (• m) c ) 

+ [l,2,l]((o) c x (o o) c x (m) c ) 

+ [2,i,i]((» •),. x (•),. x (m) c ) 

+ [l,l,l,l]((«) c x(«) c x(») c x(«) c ) 



F. Fusion rules 

While the details of each fusion rule between diagrams (see 
Eqs. (24a), (24b) and (24c) in Section II B 2) is lattice-specific, 
important and general fusion properties can be described us- 
ing the lattice-independent diagrammatic notation introduced 
in Appendix V E. Indeed, as shown in Appendix V G, these 
lattice-independent fusion rules are sufficient conditions to 
demonstrate that any cumulant given by Eq. (79) is connected, 
hence proving that H a s is local. 

When computing products, diagrams can fuse together by 
producing all possible contractions of terms Q Q between din- 
stinct blocks (••-•) x ... x (O O -•)■ Fusions inside 
blocks (• • - •) are not allowed. 

As as an example, let us apply these rules to the products 
appearing in ((• #) c ) and ((# • #) c ). 



((•) x (•)) = ((• •)) + ((•) x (•)) (84) 
= ((• •)) + ((•-•)> 



((•) x (• •)) = {(• • •)) (85) 

+ {(•jTci •)} + ((•) x (• »)) + ((iwi i)> 
= ((««»)) + 2((CKD •)} + <(*«)) 



((O O) x (o)> = ((O O •)) (86) 



+ ((• •) x (•)) + ((» •) x (•)) + ((• •) x (•)) 
= ((»»»))+2((»-» .)) + ((»-**)) 



((•) x (•) x (•)) = ((• • •)) (87) 
+ x (•)) + {{ %) x C) x (m ) 

+ ((•) x (^TS)) + ((•) x (•) x (•) 



+ ((•) X (•) X (•)) + ((•) X (•) x (•) 
+ ((•) x (6) x (•)) 

= ((O O O)) + 3((00 O)) + 3((0O0)> + ((&)) 

Combining Eqs. (84), (85), (86) and (87) with (82) leads to: 

((0) c > = ((O)) (88) 
((••),) = -((•-•)) (89) 
((0 0)J=3((»^))+2((^)) (90) 

G. Linked cluster theorem 

In this appendix we demonstrate that the result suggested 
by Eqs. (88) is general: only connected diagrams contribute 
to ((• •••••),). 

The proof will be given in two steps. First, we show that 
((O 0-0) c ) can be reexpressed as Eq. (92) by relaxing the 
constraint on internal contractions. In the second part we use 
Eq. (92) to establish Eq. (99) which shows that cumulants of 
any given order can be reexpressed with lower order cumu- 
lants hence providing a demonstration of the above mentioned 
result by direct induction. 

1. Relaxing internal contractions 

The demonstration uses the key fact that internal contrac- 
tions 

(• m^m) 

are not allowed when evaluating ((O 0-0) c ). 

Let us relax this constraint and introduce ((# •••••),,)* 
whose definition is the same as ((# •••••) ) but evaluated 
with at least one internal contraction. Obviously, the sum 

Cl m = ((••-•).) + <(••-•),)* (91) 

m m 

corresponds to a fully unrestricted evaluation where both in- 
ternal and external contractions are allowed. In that case 
each term of equation (79) produces the same set of terms by 
contraction since all partitions become equivalent when the 
distinction between internal and external contractions is sup- 
pressed. Thus, the corresponding weight can be obtained by 
the formal identification (O O - O) <H> 1 into equation (79) or, 
according to Eq. (76), z^ n ' f-> (1/nl) into equation (67). The 
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later implies Z «-> cxp(/3) and InZ «-> /3. This shows that 
fli = ((o)) and f2 m = for m > 1 which leads to, 



((O O ••••),,) 



((••-•)«.)* ifm>l, 
for m = 1. 



(92) 



<(o)> 



2. Cumulant order reduction 

The proof of the result will be given by mathematical in- 
duction. We suppose the result to be true up to rank to — 1. 
For to > 1, 

((oo-oy 

(-l) n (n- l)!x 

V [«]<(• •••••) x... x (•••«•))*. (93) 

' ' ~ ' 

Notice that the sum over n has been truncated up to m — 1 
since the n = m term contains only size-1 blocks and can 
obviously not produce any internal contraction. This equation 



m 

771— 1 

E 

71=1 



can be made explicit as, 

((••-•) x ... x (• • -••))* 

Sl 

= -<(••"••)) 



S ' 



E-E E E n„ 



9l =l Sn =l 7lGS n 7 „ S 5^P=l ffp - 11 ?= l7p - 

,/Yl 7? 7i 9 \ ,7^ 7^ 7^" 

((• •■■■•) x ... x (• m- •)) (94) 



where • denotes the result of all connected contractions be- 
tween 7 operators •. Note that the presence of the first term 
on the r.h.s of equation (94) takes into account that at least one 
internal contraction has to be performed. 

For 7 P £ Sg^ with p = 1, . . . , n and s € S™ we have 

Ep=i s P = Ep=i Ef=i l q P and thus < 

71 | 

W II ,rrt" = t^i, • • • , <?»}][7i U . . . U ln ]. (95) 

p=l 9p- \.lq=l 7p ! 

Using Eq. (95) and combining Eq. (93) with Eq. (94) leads 

to, 



<(0 0-0) c )= (96) 
v v ' 

m 

771 — 1 Sl S n 1 2 Sl 1 2 g n 

^ (-i)«(„ -i)i E-E E ■■■ E [{ffi.---.u][7iu...u 7 n]((« •)x...x(« •-•)), 

11=1 se5r _1 f 1=1 9>»= 1 7ie5g; 7neS|^ 

where summations over s have been restricted to s € to ensure that at least one internal contraction occurs. Next, a simple 

change of variables from j to fi and a sum inversion gives, 

m— 1 m— 1 

Ex — ^ x — x — Ml M2 Mgi MG 
(-!)>- 1)!^ I>M • • • < 

(97) 

m— 1 G 

E m Ec- 1 ) ( n - !) ! E • _?j x ■■•><(•••■■ •)> 

G=l fiGSp 71=1 oG5 G ' " ' 

I 



This last expression involves explicitly the cumulants defined 
in Eq. (79) and as a result, the order-™ cumulant can be ex- 
pressed as a combination of lower order cumulants: 



where the operators in the l.h.s. of Eq. (98) are arbitrary con- 
nected operators •. With T = Ep=i r P' 



<(o o-o) c ) = - V •-•)c) (98) 



G=l jugSg 



Ti r 2 t„ 
((• • « 



i> = -E Em« 

G=i ^sg 



Ml M2 MG 



)c>, (99) 



The demonstration can be readily extended to the case 



which concludes the recursion proof. 
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Interestingly, using Eq. (98), one easily recovers Eq. (88): 

<(0) c ) = ((O)) (100) 
((•©) c ) = -[2]((«) c > = -<(*0)) (101) 



((•••) c >=-[3]((») c 



(102) 



-[l,2]((i •) C >-[2,1]<(« i) c > 

((&))+ 
+ 2x^(2((w)) + ((^))) 

= 3(0 



)>+2<(^)> 



H. Fully connected diagrams 

In practice the class of fully connected diagrams, 




where each part O is connected to the others plays a very im- 
portant role when computing InZ. Indeed, they lead to the 
spatially most compact and prominent terms in Hes, as ex- 
plained in Section II B 3. The aim of this appendix is to com- 
pute the weight w m of this diagr am in the order-m. cumulant: 



((• •-•)«.) =w m ( 




(103) 



less connected terms 



Noticing that the way to connect two fully connected diagrams 
into a single fully connected one is unique, the recursion rela- 
tion for w m is readily given by Eq. (98), 

m— 1 

w ™ = - £ w p M- (104) 

p=i nes™ 
Solving this recursive relation with wi = 1 leads to, 

iy ro = (-l) m - 1 (m-l)! (105) 

Finally, the contribution to In Z, obtained by replacing 
Eq. (105) into Eq. (70) using Eq.(77), is 



InZ 




(106) 



m— 1 

less connected terms. 
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We present a full list of fusion rules (Section I) and processes (Section II) required to evaluate the effective hamiltonian 
(Section III) up to order a 14 . In Section III we also summurize the result of the resummation when it can be obtained. 

I. FUSION RULES FOR THE KAGOME LATTICE UP TO ORDER 14 

A. Order 8 

l{. Q Q } /W + 2 A>~<7 / 

B. Order 10 




II. EXPANSION OF O AND H UP TO ORDER a 1 



Processes 



O H 



Processes 



O H 



Processes 



O 



H 



Processes 



O 



u 



MM/ 



ff 

g 





.,10 



.,10 



.,10 



.,10 



.,10 



.,10 



.,10 



.,10 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



-3Ja 4 
-2Ja 6 
-2Ja 6 
-2Ja e 
-6Ja 8 
-Ja 8 
-J a 8 
-J a 8 
-5Ja 10 
-5Ja 10 
-5Ja 10 


-AJa w 
-4Ja 10 
-4Ja 10 
-4Ja 10 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 



/g 



#7 



m/ 



/A77 / 



a 12 -3Ja 12 
a 12 -3Ja 12 
a 12 -3Ja 12 
a 12 -3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 
-3Ja 12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



.,12 



/ V / 

ffi 




-3Ja 12 
-3Ja 12 
-3Ja 12 
-9Ja 12 
-4Ja 12 
-4Ja 12 
-4Ja 12 



Z 2 ^ 

/ / 



m 



-AJa 



12 



a 12 -AJa 12 
a 12 -AJa 12 
a 14 -2Ja 14 



.,14 



„14 



„14 



,14 



,14 



,14 



„14 



„14 



„14 



„14 



,14 



„14 



„14 



,14 



„14 



,44 



,44 



,44 



,44 



-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 
-2Ja 14 



''CO- 



XY/ 



/ <> <> /\, 



/<><> / / / 



^ 7 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



,44 



-2Ja 14 
-6Ja 14 
-6Ja 14 
-6Ja 14 
-6Ja 14 
-6Ja 14 
-5Ja 14 
-5Ja 14 
-8Ja 14 
-8Ja 14 
-8Ja 14 
-3Ja 14 
-7Ja 14 
-7Ja 14 
-7Ja 14 
-7Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
-3Ja 14 
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III. EXPANSION OF Km UP TO ORDER a 14 AND FULLY RESUMMED SERIES 



The symbol means that the resummation has not been performed. 

It is interesting to note that for a given dimer configuration and for a given hexagon or pair of neighboring hexagons of the 
lattice there is one and only one active kinetic process in the list. 



Processes 



Order 14 



Processes 



Order 14 



DO 



Processes 



Order 14 



Processes 



Order 14 



DO 



/ST / 



si 



/ST / 

& 



xs, 



M 
M 
M 

/&/ 



-3 Ja 4 (a 8 + 1) 
3Ja 8 
-2Ja 6 
-2Ja 6 
-2Ja 6 
2Ja 12 
2Ja 12 
2Ja 12 
-Ja 8 
-Ja 8 
-Ja 8 











3Ja 10 



3Ja 



10 



2 

3Ja 10 



7Ja 



7Ja 



7 J a 



ll 



3.7a 4 



3Ja B 



2,7a D 



2Ja° 
a 12 -l 

2Ja 6 
a 12 -l 

2Ja 



12 



Q 12 -l 
n 7 12 

_ 2Ja 
a 12 -l 

n 7 12 

_ 2Ja 
a 12 -l 

Ja 8 
a 16 -l 

Ja 8 

Q 16 -l 

Ja 8 
a 16 -l 

7 16 
a 16_! 

7 16 

Ja 16 
ai 6 -l 





V 





*7 7 14 

7Ja 

8 

-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 






AAV 



-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 
-Ja 12 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 




Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



,7a 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
2 



Ja 1 
4 



,7a 1 
4 



Ja 1 
4 



Ja 1 
4 



Ja 1 
4 



3JU 1 



3.7a 1 
4 



